Quantum limited particle sensing in optical tweezers 
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Particle sensing in optical tweezers systems provides information on the position, velocity and force of the 
specimen particles. The conventional quadrant detection scheme is applied ubiquitously in optical tweezers 
experiments to quantify these parameters. In this paper we show that quadrant detection is non-optimal for 
particle sensing in optical tweezers and propose an alternative optimal particle sensing scheme based on spatial 
homodyne detection. A formalism for particle sensing in terms of transverse spatial modes is developed and nu- 
merical simulations of the efficacy of both quadrant and spatial homodyne detection are shown. We demonstrate 
that an order of magnitude improvement in particle sensing sensitivity can be achieved using spatial homodyne 
over quadrant detection. 



I. INTRODUCTION 

The application of the radiation pressure force for the trap- 
ping of atoms and neutral particles was pioneered by Arthur 
Ashkin JD]. This was followed by a plethora of seminal ex- 
periments utilizing the radiation pressure force 101, for ex- 
ample in the displacement and levitation in air and water of 
micron-sized particles and together with Steve Chu, for 
the development of a stable three-dimensional atom cooling 
and trapping experiment using frequency-detuned counter- 
propagating laser beams fl]. In particular, the demonstration 
of optical tweezers JH], based largely on the transverse gra- 
dient force of a single focused Gaussian optical beam was a 
significant contribution to optical trapping in biology [6]. 

In biological systems, optical tweezers were first used to 
trap and manipulate viruses and bacteria 0]. This was fol- 
lowed by a burgeoning number of experiments using optical 
tweezers for measurements of DNA/RNA stretching and un- 
folding i^lHd^Gjldsl], intracellular probing, manipulation 
of gamete cells, trapping of vesicles, membranes and colloids 
[Qj, Q and DNA sequencing using RNA polymerase lfl5ll . 
In particular, for the first time, quantitative biophysical stud- 
ies of the kinetics of molecular motors ifToTl (e.g. myosin lfl7tl 
and kinesin [18]) at the single molecule level was made pos- 
sible with the use of optical tweezers. Coupled with conven- 
tional position sensitive detectors (i.e. using quadrant pho- 
todetectors lfl3i [T9i l2(lo . the position of, and force on, a 
bead tethered to a molecular motor can be measured at the 
single molecule level l2ll I22I l23ll . The sensitivities attain- 
able for force and position measurements of particles in op- 
tical tweezers are in the sub-piconewton and sub-nanometer 
regimes Il3l[l7il23ll . respectively. The application of the op- 
tical tweezers technology has led to a more complete bio- 
physical understan ding of the kinetics of molecular motors 

in m m m m mr m eh m in m m - a quintessen- 
tial demonstration of new physical techniques yielding new 
insights into biology. 

Beam position and momentum sensing is particularly cru- 
cial for particle sensing in optical tweezers enabling high- 
precision particle position and force measurements jl3l [jjj 
12011 . Therefore it is important that such measurements are 
performed optimally to achieve the highest measurement effi- 



cacy. Recently, Hsu et al. [[36[] showed that the conventional 
quadrant detection scheme is non-optimal for measurements 
of the position and momentum of optical beams, even in the 
absence of classical noise sources. An alternative scheme for 
the optimal detection of the position and momentum of an op- 
tical beam was proposed, based on a spatial homodyne detec- 
tion scheme. This scheme has also been proven to perform at 
the quantum limit of light based on Cramer-Rao informational 
bounds rf37ll . Therefore, it has become apparent that the use 
of quadrant detection for particle sensing in optical tweezers 
systems is non-optimal; and the introduction of spatial homo- 
dyne detection could offer the possibility for greater particle 
tracking sensitivities. 

In this paper, we address the pertinent questions for parti- 
cle sensing technology in optical tweezers systems - have we 
reached the limit of particle tracking sensitivity and can this 
limit be surpassed using quantum resources? We believe that 
in answering this technique related question, naturally arises a 
biophysical question - i.e. with significantly enhanced sensi- 
tivities, are we able to detect molecular kinetics, at the single 
molecule level, that were previously unresolvable? This bio- 
physical question has wide implications as there are many vi- 
tal protein conformational changes that occur in the angstrom 
regime, and within millisecond timescales ll23"ll . For exam- 
ple, molecular motors move along nucleic acids in steps of a 
single-base pair scale (e.g. 3.4 A on dsDNA) l23ll and the bac- 
terial DNA translocase FtsK moves at speeds of 5 kilobases 
per second (HI. Therefore, enhanced particle sensing could 
elucidate these finer features with greater sensitivity than con- 
ventional particle sensing techniques in optical tweezers sys- 
tems. 



This paper begins by formalizing an optimal parameter esti- 
mation procedure for particle sensing based on the analysis of 
the spatial properties of the field scattered by a particle in an 
optical tweezers. We show that split detection is non-optimal 
and consequently propose an optimal measurement scheme 
based on spatial homodyne detection. The efficacy of parti- 
cle sensing is evaluated using the signal-to-noise ratio (SNR) 
and sensitivity measures; and the efficacy of spatial homodyne 
detection and split detection systems are compared. 
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II. OPTIMAL PARAMETER ESTIMATION FOR SPATIAL 
PROPERTIES OF OPTICAL FIELDS 



a sum of coherent amplitude components and quantum noise 
operators, given by 



An optical field can be formalized and described using a 
range of parameters - e.g. the polarization, the amplitude- 
phase quadratures, and the transverse spatial profile. These 
parameters can be measured using a range of detection tech- 
niques (e.g. polarimetry, direct detection, interferometry and 
beam profiling) and an estimate of their values in the presence 
of classical and quantum noise, and detection inefficiency is 
obtained. Here we develop a formalism to quantify an arbi- 
trary spatial modification of the field parameterized by a pa- 
rameter p (e.g. p could quantify the displacement of a spa- 
tial mode along a transverse axis l36lo . In principle, an arbi- 
trary field can be treated and the field properties can be mod- 
eled using Maxwell's equations ll39ll . However, for spheri- 
cal fields such as those produced by scattering processes from 
small particles, after optical imaging of the field, the paraxial 
approximation is valid and the propagating field can be de- 
scribed using two-dimensional spatial modes in a convenient 
basis. 

The sensitivity of measurements on optical fields is ulti- 
mately limited by quantum noise on the fields, exhibited typ- 
ically as shot noise. To understand such limits it is impor- 
tant to use a full quantum mechanical description of the field. 
The spatial quantum states of an optical field exist within an 
infinite dimensional Hilbert space. Depending on the spa- 
tial symmetry of an imaged optical field, the spatial states 
of the field may be conveniently expanded in the basis of 
the rectangularly-symmetric TEM.„ m or circularly-symmetric 
LG„.z modes. A field of frequency uj can be represented 
by the positive frequency part of the electric field operator 
E + (p)e Mt . We are interested in the transverse information 
of the field described fully by the slowly varying field enve- 
lope operator E + (p), given by 



E+(p) = l 



h,uj 



>(p), 



(i) 



where p = (x, y) is a co-ordinate in the transverse plane of 
the field, and the summation over the parameters j, m and n 
is given by 



E = E EE- 

,7,771,71 j£{x,y} m — o n — o 



(2) 



In this paper, we adopt the TEM,„„ mode basis for con- 
venience, such that u J mn (p) and a? mn are, respectively, the 
transverse beam amplitude function and the photon annihila- 
tion operator for the TEM mn mode with polarization j. The 
u m „(p) mode functions are normalized such that their self- 
overlap integrals are unity, so that the inner product 



(P) (p) 



,{p)dp 
(3) 



We now apply an arbitrary spatial perturbation, described 
by parameter p, to the field. Eq. ([TJ can then be rewritten as 



E+(p,p) 



flUJ 
TlUJ 

2^V 



J, 771, 71. 



(4) 



where 



a(p)v(p,p) = ^ {aP mn )-a> mn {p,Q) = ^ {aP mn )u mn {p, 0)j, 



a(p) is the coherent amplitude of mode v(p,p), and j is the 
unit polarization vector. We see from Eq. (HJi that a(p) and 
v(p, p) can be related to ~E + (p,p) by 



a(p) 
v(p,p) 



'^{p,p),^{p,p)) 



flLO 



-iN v B (p,p) 



(6) 



where E + (p,p) = (E + (p,p)), and the normalization con- 
stant N v is given by 



E + (p,p),E + (p,p) 



-1/2 



E" r (p,p) -E + (p,p)dp 



-1/2 



(7) 



Note that |a(p)| 2 is the mean number of photons passing 
through the transverse plane of the field per second and in 
this paper we assume a(p) to be real, without loss of gen- 
erality. The quantum noise operator corresponding to mode 
(P) (p,0) is given by 8a> mn . 

In the limit of small estimate parameter p, the Taylor ex- 
pansion of the first bracketed term in Eq. (|4|i is given by 



a(p)v(p,p) w a(0)v(p,0) +p 



d[a(p)v(p,p)] 



dp 



(8) 



p=0 



where the first term on the right-hand side of Eq. dS) indicates 
that the majority of the power of the field is in the v(p, 0) 
mode. The second term on the right-hand side of Eq. ([8]) de- 
fines the spatial mode w(p) corresponding with small changes 
in the parameter p 



w(p) 



d[a(p)v(p,p)} 



dp 



(9) 



p=0 



From Eq. ^ we see that the amplitude of mode w(p) is di- 
rectly proportional to the magnitude of the spatial perturbation 
of the field. 
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A. Split detection 



B. Spatial homodyne detection 



For optical beam position and momentum measurements, 
the conventional detection scheme used is split detection (a 
one-dimensional quadrant detector). In split detection, the op- 
tical beam under interrogation is incident centrally on a split 
detector, as shown in Fig.Q](c). The difference between the 
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FIG. 1 : (a) Schematic diagram and (b) wavefront illustration for an 
optical tweezers experiment. A trapping field is focused onto a parti- 
cle. The particle scatters the incident trapping field, with the resulting 
scattered and residual trapping fields collected by an objective (Obj) 
lens. This is followed by imaging of the collected fields onto a posi- 
tion sensitive detector (PSD) in the far-field using an imaging (Img) 
lens. The position sensitive detector could consist of either a (c) split 
detection system or a (d) spatial homodyne scheme with the local 
oscillator (LO) beam in an optimized spatial mode for the relevant 
measurement of parameter p. 



photocurrents from the two halves of the split detector con- 
tains partial information about the position/momentum of the 
beam, given by l36ll 



Ai 



SD 



2epV 
Tiuj 



/0 poo 
E +t ■ E + dp - E +t ■ E+dp 
-oo JO 



(10) 



where X^ = a j, + a/ is the amplitude quadrature operator of 
the flipped mode with transverse mode amplitude function 

-v(p,0), x<0 (U) 

The amplitude quadrature operator can be written in terms of 
its coherent amplitude 



a f (p) = a(p) (v/(p),v(p,p)) 



(12) 



wherein resides the signal due to the parameter p, and a quan- 
tum noise operator SX^ = X^ — (Aj ) which is ultimately 



responsible for placing a quantum limit on the measurement 
sensitivity so that 



Ai SB = a(p) 2a f (p) 



■SXj 



(13) 



The (v/(p), v(p,p)) term in Eq. (fT2l is the overlap inte- 
gral between the flipped mode v/ (p) and the displaced mode 
v(p,p). 



Hsu et al. [ 36] proposed a new displacement measurement 
scheme that is optimal for detecting beam position and mo- 
mentum. The spatial homodyne scheme utilizes a homodyne 
detection setup that has a local oscillator mode optimized for 
the displacement measurement of the input beam, as shown 
in Fig. [TJ(d). The local oscillator (LO) beam interferes with 
the input beam on a 50/50 beam-splitter. The outputs of the 
beam-splitter are then detected using a pair of balanced single- 
element photodetectors, with the difference in photocurrents 
providing the measurement signal. The spatial homodyne 
scheme was also proven to perform at the Cramer-Rao bound 
rf37ll . therefore extending the capabilities of the spatial homo- 
dyne scheme for the optimal measurement of any spatial pa- 
rameter p (e.g. the measurement of the orbital angular mo- 
mentum of light Hot]). We now proceed to derive the pho- 
tocurrent for the spatial homodyne detection scheme. 

The input beam (as described in Eq. (0]i) is interfered with 
the bright LO beam with mode-shape w(p). The positive fre- 
quency part of the electric field operator for the LO given by 



eJoCp) 



"LOw(p) 



E 

J , 771 , 71 



(p) 



,i(j> 



(14) 



where is the phase difference between the local oscillator 
and the input beam. 

The photocurrent at each photodetector (distinguished by 
the subscripts + and -, respectively), assuming detectors of 
infinite extent, is given by 



2e V 
2e V 



(15) 



(E+±E+ )t.(E+±E+ )dp, 



(16) 



whereby one output of the spatial homodyne attains a it phase 
shift with respect to the other output due to the hard-reflection 
from the beam-splitter. 

Substituting Eqs. ([TJ) and ( fT4b into Eq. ( [Tol l and taking the 
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subtraction of the photocurrent from the two detectors gives 

POO 



= CCLO 



+e l0 w(p)- 5 L U L(P.P) 

Smn( w (p) 



r/p 



3,m,n 



(17) 



where o w is an annihilation operator describing the compo- 
nent of the input field in mode w(p), and by definition the 
= e~ l ^d w + e^alu is the quadrature operator of that 
component at phase angle <\>. In the above, we have taken the 
condition a^o ^> (a w ) and invoked the linearization approxi- 
mation, thereby removing terms that do not involve oilo The 
orthonormality property of modes given in Eq. (0 has also 
been used. 

An optimal estimate of the parameter p is obtained when 
the local oscillator mode w(p) is chosen to match the asso- 
ciated input mode v(p,p), as shown in Eq. ill) . The spatial 
homodyne detection scheme then extracts from the signal field 
a quadrature variable associated with the local oscillator field 
mode, with quadrature phase angle given by 0. 

It should be noted that Delaubert et al. 113711 have shown that 
optimal parameter estimation can be achieved using a pho- 
todetector array for the cases where the signal field E + (p, p) 
is shot noise limited or single mode squeezed, so long as the 
array resolution is sufficiently small. Array detection is re- 
stricted to amplitude quadrature detection and is not polariza- 
tion resolving, however in situations where these restrictions 
are satisfied is formally identical to spatial homodyne detec- 
tion. 



C. Quantifying the efficacy of parameter estimation 

We now introduce the SNR and sensitivity measures for the 
spatial homodyne and split detection schemes. 

For the spatial homodyne detection scheme, the measured 
signal is the mean signal component of the difference pho- 
tocurrent in Eq. (fTTt , given by 



(A lSH ) = a w a w (p) (e tM + e" lA *) 



(18) 



where a w (p) = a(p) (w(p), v(p,p)). For matched local os- 
cillator and signal phases such that = 0, the maximal signal 
is obtained, given by 



(At 



sh; 



2a-Loa w (p). 



(19) 



The corresponding noise component is given by 



(A; 



SH/ 



(A* SH > 2 = a LO AXt 



(20) 



where A 2 X^ = ((SX^) 2 ) is the variance of the signal field 
mode. The resulting SNR is given by 



SNR S H 



2a w {p) 
AXt° 



(21) 



If the optical field is in a coherent state, as is typical of a low 
noise laser AX^° = 1 and the SNR for the spatial homo- 
dyne detection scheme is given by 



SNR 



SH,coh 



2a w (p). 



(22) 



Clearly, although experimentally challenging, squeezing the 
signal mode such that AX^ LO < 1 has the capacity to further 
enhance the SNR. 

Alternatively, we introduce the sensitivity S measure, 
which is defined as a change to parameter p required to pro- 
vide SNR = 1 for a signal field in a coherent state, given by 



SH.coh 



9SNR 




-l 

1 


da w (p) 




dp 


p=0 


~ 2 


dp 


p=0 



(23) 



For comparison, the corresponding SNR for the split detec- 
tion scheme in the coherent state limit is given by 



SNR SDiCoh = 2a f (p), 
with a sensitivity given by 



<5sD,coh — ^ 



datfip) 




dp 


p=0 



-I -1 



(24) 



(25) 



III. PARTICLE SENSING IN OPTICAL TWEEZERS 

Fig.[TJ(a) shows a typical optical tweezers setup. A trapping 
beam in the TEMoo mode is focused onto a scattering particle. 
In this instance we assume that the particle is spherical, with 
a permittivity greater than that of the medium, €2 > ei- If 
the particle has a diameter larger than the wavelength of the 
trapping beam, light rays are refracted as they pass through 
the particle, as shown in Fig. 2. This refracted light results in 
an equal and opposite change of momentum imparted on the 
particle. Due to the intensity profile of the beam, the outer ray 
is less intense than the inner ray. Consequently, the resulting 
force acts to return the particle to the center of the trapping 
beam focus 

The effective restoring/trapping force is due to two force 
components - (i) the gradient force F grac i resulting from the 
intensity gradient of the TEMoo trapping beam, that acts trans- 
versely toward the high intensity region and (ii) the scatter- 
ing force F scat resulting from the forward-direction radiation 
pressure of the trapping beam incident on the particle. In the 
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FIG. 2: Illustration showing a TEM U0 trapping beam impinging on a 
spherical scattering particle. Rays 1 and 2 are refracted in the spher- 
ical particle, thereby undergoing a change in momentum. A corre- 
sponding equal and opposite change in momentum is imparted on 
the particle resulting in the particle being attracted to the center of 
the trapping beam. F grat j and F scat are the gradient and scattering 
forces, respectively, ei and are the respective permittivity of the 
medium and the sample. 



B. Propagation of fields through system 

We now formalize all the relevant fields that propa- 
gate through the optical tweezers system, as shown in the 
schematic of the optical tweezers arrangement of Fig. Q](a). 
Fig.[T](b) illustrates the wave-front of the trapping and scat- 
tered fields. The trapping field is incident from the left of the 
diagram and is then focused onto a spot, from the focusing 
lens. The particle is trapped near the center of this focal spot 
and scatters the incident trapping field, with the forward scat- 
tered and residual trapping field being collected by the objec- 
tive lens. This is followed by imaging into the far-field onto a 
position sensitive detector. 



Trapping field 



focal region of the optical tweezers trap the gradient force is 
typically dominant. 

It is important to note that in some optical tweezers experi- 
ments the trapped particle has radius less than the wavelength 
of the trapping laser. In this regime, the trapping force on the 
particle is generated due to an induced dipole moment. The 
dipole moment induced will be along the direction of trapping 
beam polarization. The assumption that the particle is spher- 
ical is no longer important, since the particle has no struc- 
tural deviations greater than the wavelength of the trapping 
beam. This allows the particle to be treated as a normal dipole, 
hence the particle experiences a force due to interaction of its 
induced dipole moment with the transverse electromagnetic 
fields of the impinging light. This force is proportional to the 
intensity of the beam and has the same net result as before; it 
acts to return the particle to the center of the trapping beam 
focus. 

The position and force sensing of the trapped particle can 
then be obtained by imaging the scattered field from the par- 
ticle on a position sensitive detector such as the commonly 
utilized quadrant photo-detector |[T3i[T9ll20ll . or a spatial ho- 
modyne detector !^ . 



A. System configuration 

The collection efficiency of the light field is given by the 
numerical aperture (NA) of the objective lens (as shown in 
Fig.Q](b)), given by 



Assuming that the trapping field is Gaussian and hence in a 
TEMoo mode, the positive frequency part of the electric field 
for the trapping beam at the waist of the trap (denoted by the 
superscript T), is given by 



2^V 



a trap uj (p), 



with the mode-shape function given by 



u oo(p) 



Wry/ft 



2 / 2 

3 9 ^Ptrap, 



(28) 



(29) 



where p 2 = |p| 2 , Wt is the waist size of the trapping beam, 
and p trap is a unit vector representing the polarization of the 
trapping field. 

Using the paraxial approximation, the positive frequency 
part of the electric field of the trapping beam after propagation 
of a distance /o from the focus to the objective lens (of focal 
length fo) is given by 



E trap(P) = * 



flUJ 



^a trap e-^°u°n fl (p), (30) 
ZcqV 



where k = 2it/ X is the wave-vector of the trapping field. With 
the exception of the replacement wo — > wt, Uq (p) is defined 
identically to uj , with the radius of the spot at the objective 
being wq given by 



w 



/oA 
ttwt 



(31) 



NA = nsm.6. 



(26) 



where n and 9 are the refractive index and the collection half- 
angle of the lens, respectively. 9 is related to the lens diameter 
D (assuming the object is at the focus, with focal length f£) 
by 



tan 9 = 



D/2 



(27) 



Aperturing due to the finite radius R of the objective lens is 
taken into account via the aperture function IIr(p) given by 



n*(p) = 



1, p < R 
0, P>0, 



(32) 



where R can be related to the numerical aperture (NA) of the 
imaging system and refractive index of the trapping medium 

nby R = f NA/\/n 2 - NA 2 . 
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Scattered field 



Detection 



In principle, there could be multiple inhomogeneous par- 
ticles within the optical tweezers focus, scattering the input 
trapping field. For this scenario, several numerical methods 
exist to calculate the scattered field - e.g. the finite differ- 
ence frequency domain and T-matrix hybrid method fiHl and 
the discrete-dipole approximation and point matching method 
ll42tl . However, for simplicity we consider the scattering from 
a single spherical, homogeneous particle with diameter much 
smaller than the wavelength. The resulting scattered field can 
be modeled as dipole radiation, having a positive frequency 
electric field 1I39I1 given by 



Since the total field after the objective lens consists of both 
the scattered field and the residual trapping field, we now in- 
clude both fields to describe the total field after the objective 
lens, given by 



Etotal(P.P) = E scatO>P) + E trap(p) 



(40) 



After the objective lens, the beam is focused onto a detector 
in the far-field image plane, via the use of an imaging lens. 
Assuming the lens is thin and ideal, the field in the image 
plane is obtained by taking the Fourier transform of Eq. d40l i, 
given by 



E Lt( r ,P) 



ei - £2 
ei + 2e 2 



-r' x f x E trap (p), 



r 



E 



total 



(33) 

where r = lie + yy + zi is the coordinate of the field with 
respect to the center of the optical tweezers, r' = (x — p)x + 
yy + zi is the coordinate of the field with respect to the dis- 
placed particle, r' = |r'|, and f = r'/r'. The radius of the 
spherical scattering particle is given by a. 

The scattered field is then collected by the objective lens 
(as shown in Fig.Q](b)), with the corresponding positive fre- 
quency part of the electric field given by 



^( E taIl(P.P) 
■0+ 



(41) 



= .FIE 

= E^ t (r, P ) 



S cat(P,P)) + ^(E trap (p)) (42) 



E 1i,(r) 



(43) 



E°cat(P,P) = 



'E + 



scat 

(r ,p) • 1 1 + E 

scat 

(r ) '&n 




(34) 



where T = (X, Y) are the transverse co-ordinates in the im- 
age plane. It is important to note that the analysis presented 
here is independent of the absolute scaling of the image plane 
co-ordinates. In an experimental situation a scaling factor is 
introduced that depends on the choice of magnification lenses 
used. 

The critical parameters for assessing sensitivity of parti- 
cle monitoring are a(p), v(T,p) and w(T). These parame- 
ters can now be calculated using Eqs. (0, © and (|9j. Using 
Eq. (0 we now find 



eoc y r Q , 

(to' x rb' x u£ ) i + (to' x tq' x u£ ) n 



a(p) 



(35) 



/ E tl f ap (r),E t J f ap (r) 



CCtrap , 



(44) 

(45) 
(46) 



where ro = xx + yy + foi and r' 
The unit vectors 1, m, and n 



■p)x + yy + /oz- 



1 



-f(y,p-x,o) 
t> 

— — {-fo(x -p),-foy,p' 2 ) 

pro 

—r(x-p,y,0), 
P 



(36) 
(37) 
(38) 



are used to include the effect of the objective lens on the polar- 
ization of the scattered field, where p' = (x — p)x + yy and 
p' = \p'\. The term -J fo/r' Q describes the compression of 
the intensity of the scattered field due to the change in propa- 
gation direction induced by the objective lens. To simplify the 
equation, we have defined the constant K given by 



where we have assumed that the trap power is greater than 
the scattered power, as is the case for scattering from a small 
particle; and for simplicity that only the scattered field is aper- 
tured by the objective lens. The latter assumption is reason- 
able for optical tweezers systems with a sufficiently large trap 
waist size and numerical aperture. In this paper, we restrict 
our analysis to the realistic scenario of NA > 0.2 and choose 
a trapping field waist size of 4 /im. With these parameters, 
trap field clipping due to the aperture causes only 15 ppm loss 
and is therefore negligible. 
Using Eq. (|6]l we obtain 



v(T,p) 



ap 



■ E taap( r ~ 



(47) 



2„3 



K = CKtrapfc a 



ei - 62 
ei + 2e 2 



(39) 



where we have used the relations for N v and a(p) given in 
Eqs. (IT) and (|46*| |. respectively. 
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Now using Eq. (O we obtain the functional form for the 
mode that contains information about the particle position, 
given by 



w(T) 



dp 



(48) 



p=0 



Note that this mode is only dependent on the scattered field. 

We now calculate the SNR of the spatial homodyne and 
split detection schemes for particle sensing in an optical 
tweezers arrangement. Substituting the expressions obtained 
in Eqs. (gSJ) - (|48]) into Eq. the SNR for the spatial ho- 
modyne detection scheme is given by 

SNR SH ,coh = 2& w (p) 

= 2a(p)(w(r,p),v(r,p)) 



-2i 



2e V 

TlUJ 



w(T, P y 



E 



■i+ 



= -2i 



:at(i» + E trap (r)Jdr 

/oo 
w(r)*.E s + t (r, P )dr 
-oo 



2eoV 



/oo 
w(r)* ■ A(r)dT, (49) 
-oo 



where the effective aperture function in the image plane co- 
ordinates is given by 




fo e _ifc ( ro_r o-/o) 



A(r) 



In a similar manner using Eq. (I24I) . the SNR for the split 
detection scheme is given by 



(r~o' x rb' x u£ ) n 



r x r x Uq ) 1 

n B (p)Y (50) 



SNR SD 



coh 



/oo 
v/(r)* • A(r)dT (51) 
-oo 



Correspondingly, the sensitivities for the spatial homodyne 
and split detection schemes can be conveniently calculated us- 
ing Eqs. d23l l and d25l l. respectively. 



IV. SIMULATION AND RESULTS 

A formal description for the trapping and scattered fields in 
an optical tweezers configuration was presented in SectionlllTl 
We now numerically solve for the scattered field from a parti- 
cle trapped in the optical tweezers. We utilize the field imag- 
ing system shown in Fig. Q] (b) to image the scattered field 
into a propagating optical beam that is subsequently detected. 
We compare the SNR and sensitivity of both split and spatial 
homodyne detection schemes (described in SectionllTli. 

As mentioned in the preceding section, the origin of the co- 
ordinate system is defined to be at the focal point of the optical 
tweezers focusing lens system. The optical fields propagate in 



the z direction and the scattering particle was assumed to be 
spherical and homogeneous. We model particle displacement 
in the x-y plane, to illustrate the effect on the scattered field 
in the transverse plane. The far-field intensity distribution ar- 
riving at the detector is given by the interference between the 
trapping and forward scattered fields calculated from Eq. (|43l l 
and shown in Fig. [3] As the trapping field is far more in- 
tense than the scattered field, we have subtracted its intensity 
from the images shown in this figure as well as subsequent 
figures, to make visible the interference fringes between scat- 
tered and trapping fields. Note that the terms due to just the 
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FIG. 3: Interference pattern of the trapping and forward scattered 
field in the far-field image plane for 200 mW trapping power, A = 
1064 nm, particle radius a = 0.1 ^im, permittivity of the medium 
ei = 1, permittivity of the particle £2 = 3.8, and objectives with NA 
= 0.99 and focal spot size of 4 (im. We assume absorptive losses in 
the sample are negligible. Figures (a)-(c) and (d)-(e) assume the trap- 
ping field is linearly x and y-polarized, respectively. The color bar 
shows scale of the intensity distribution. The particle displacements 
are given by (a), (d): 1 jj,m; (b), (e): 0.5 [im; and (c), (f): /im. 



scattered field have been ignored to reduce numerical error, 
justified since the total scattered power is four orders of mag- 
nitude smaller than the trapping beam power. The detection 
area was chosen to be larger than the area of the calculated 
image field to avoid inaccuracies due to clipping of the image. 
Notice that as the particle moves in one direction, the intensity 
distribution shifts in the opposite direction, due to the lensing 
effect of the objective. Note also the difference in intensity 
distribution between the x and y trapping beam polarization 
directions - i.e. the interference pattern appears "compressed" 
along the polarization axis due to the dipole scattering distri- 
bution of the particle. 

The SNR for the split and spatial homodyne detection 
schemes were calculated, the results of which are shown in 
Figs. [4] (a)-(c). The SNR of the split detection scheme was 
evaluated by applying Eq. ( BTT l to the calculated interference 
signal, shown in Fig.|4](a). To calculate the SNR for the spa- 
tial homodyne detection, the optimal LO mode first had to be 
determined. Improved SNR is possible with spatial homodyne 
detection when compared with split detection for all particle 
displacement regimes. However, the optimal LO mode de- 
pends on the position of the particle, so to achieve this a dy- 
namical mode optimization routine would need to be imple- 
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-1 -0.8 -0.6 -0.4 -0.2 0.2 0.4 0.6 0.8 1 
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FIG. 4: Normalized SNR versus particle displacement for (a) split 
detection, (b) spatial homodyne detection with LO spatial mode op- 
timized for small displacement measurements and (c) spatial homo- 
dyne detection with LO spatial mode optimized for larger displace- 
ment measurements. The black solid and red dashed lines are for 
linearly x and y-polarized trapping fields, respectively. The LO spa- 
tial modes for the small displacement measurements are (d): y and 
(e) a>polarized trapping fields, whilst for large displacement mea- 
surements are (f): y and (g) a>polarized trapping fields. 



merited. Here we present results with detection optimized for 
two specific cases: (i) for a particle located close to the ori- 
gin (p <§C wt) as modeled in the theory section; and (ii) for a 
particle displaced from the origin by a factor of order wt- 

For the small displacement limit, the LO field was deter- 
mined from the first order term in the Taylor expansion of 



Eq. (0 for the scattered field. The resulting SNR is shown 
in Fig. |4] (b); with the corresponding LO spatial modes as- 
suming y and x linearly polarized trapping fields shown in 
Figs. [4] (d) and (e), respectively. One observes that for dis- 
placements significantly less than the trapping beam waist size 
the SNR is linear, with the optimum sensitivity - correspond- 
ing to maximum slope in the SNR - occurring at zero displace- 
ment and significantly surpassing that achievable with split 
detection. Particle tracking with optimum sensitivity is possi- 
ble in this linear regime. At particle displacements of around 
~ 1 0.4 1 /xm, however, the SNR peaks. Small displacements 
of a particle around these points leave the SNR unchanged. 
Hence the signal read out from the spatial homodyne detector 
also remains unchanged, with the result that particle tracking 
becomes ineffective. As the particle position increases fur- 
ther, it moves out of the trapping field, causing a drop in the 
total scattered power and consequential exponential decay in 
the SNR. 

It is possible to recalculate the LO field mode to optimize 
the sensitivity for particles fluctuating around any arbitrary 
position by performing a Taylor expansion in p of the scat- 
tered field about that position, and retaining only the first order 
term. Fig.|4](c) shows the resulting SNR and corresponding 
LO mode shapes when the LO mode is optimized for parti- 
cles fluctuating around 0.4 /im. Notice that now the maxi- 
mum SNR slope, and hence optimum sensitivity, is shifted 
from zero displacement to displacements of around 0.4 [im. 
Hence, we see that as the tracked particle moves, it is possi- 
ble to dynamically adjust the LO field shape to optimize the 
measurement sensitivity and hence the particle tracking. 

We now numerically evaluate the sensitivities of the split 
and spatial homodyne schemes in the small displacement 
limit, given by Eqs. d25l l and (l23l respectively. The sensi- 
tivity is the minimum detectable displacement, defined as the 
displacement required to change the SNR by 1. The respec- 
tive sensitivity curves for (i) split and (ii) spatial homodyne 
detection versus the numerical aperture of the objective lens 
are shown in Fig. [5] 

The minimum detectable displacement for both the split 
and homodyne detection schemes decrease with increasing 
NA of the collection lens. As the NA increases, more of 
the scattered field is collected, therefore providing more in- 
formation about the scattering particle. The spatial homo- 
dyne outperforms the split detection scheme for all NA val- 
ues. This is due to the spatial homodyne scheme providing 
optimal information extraction of the detected field whereas 
the split detection scheme only measures partial information 
of the detected field, as derived in Eq. ( fT2] i. Therefore curve 
(ii) is the quantum limit for particle sensing in optical tweez- 
ers systems. In order to perform measurements below this 
quantum limit, non-classical resources have to be used. For 
example, squeezed light in the spatial mode ll43ll correspond- 
ing to the displacement signal mode can be injected into the 
optical tweezers system to reduce the quantum noise floor and 
therefore enhance position sensing (36]. 
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FIG. 5: Mininmum detectable displacement versus collection lens 
NA for (i) split and (ii) spatial homodyne detection, normalized by 
K. The solid and dashed lines are for linearly x- and y-polarized 
trapping fields, respectively. The axis on the right shows the min- 
imum detectable displacement assuming 200 mW trapping power, 
A = 1064 nm, particle radius a — 0.1 /im, permittivity of the 
medium ei = 1, permittivity of the particle £2 = 3.8, and objec- 
tives with focal spot size of 4 /im. We assume absorptive losses in 
the sample are negligible. The split detection non-optimality shaded 
area shows the particle sensing sensitivity loss due to incomplete 
information detection from split detection. The quantum resources 
shaded area indicates the region where quantum resources such as 
squeezed light l43ll can be used to further enhance the sensitivity of 
particle sensing measurements. 



V. CONCLUSION 



We have developed a formalism for particle sensing in op- 
tical tweezers via the analysis of the transverse spatial modes 
imaged from a scattered field. The conventional quadrant de- 
tection scheme, used ubiquitously in optical tweezers experi- 
ments, was shown to only detect partial information from the 
scattered light field. We propose instead the use of spatial ho- 
modyne detection whereby optimal information from the scat- 
tering particle can be obtained via the appropriate transverse 
spatial mode-shaping of the LO field. A numerical simulation 
of the SNR and sensitivity of both split and spatial homodyne 
detection was presented and we demonstrate that up to an or- 
der of magnitude improvement in the sensitivity of spatial ho- 
modyne over split detection can be achieved. 
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